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Abstract 

Branching process approximation to the initial stages of an epi- 
demic process has been used since the 1950's as a technique for pro- 
viding stochastic counterparts to deterministic epidemic threshold the- 
orems. One way of describing the approximation is to construct both 
branching and epidemic processes on the same probability space, in 
such a way that their paths coincide for as long as possible. In this pa- 
per, it is shown, in the context of a Markovian model of parasitic infec- 
tion, that coincidence can be achieved with asymptotically high prob- 
ability until M/v infections have occurred, as long as Mn = o(iV 2 / 3 ), 
where N denotes the total number of hosts. 



1 Introduction 

The classical law of large numbers and central limit theorem have process 
analogues for many Markovian models arising in population ecology. The 
law of large numbers is replaced by a deterministic process, obtained by 
solving an appropriate system of ordinary or partial differential equations, 
and the central limit theorem is replaced by a diffusion approximation around 
the deterministic limit. For many techniques and examples concerning such 
density dependent Markov population processes, see Kurtz (1976, 1981). 
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In the context of invasion biology, when the central question is whether 
the introduction of a small number of individuals of a species can lead to its 
becoming established in a new habitat, these large population approxima- 
tions are no longer appropriate. The more natural process approximations, 
at least if spatial restrictions on mixing are not critical in such small popu- 
lations, are now branching processes. These were introduced, in the context 
of epidemic theory, by Whittle (1955), Kendall (1956) and Bartlett (1956, 
p. 129); here, infected individuals play the part of the invading species, and 
those that are infected by an individual correspond to an individual's 'off- 
spring'. 

When considering the development of a single species as a branching 
process, the biological quantity Ro, the lifetime mean number of offspring 
of a single individual when unhampered by competition from others of the 
same species, is just the mean offspring number of the corresponding Galton- 
Watson process. The branching process criticality theorem then corresponds 
to the biological meta-theorem, that an invading population can only be- 
come established if its R (in the context that it experiences upon invasion) 
exceeds 1. For models involving more species, the analogy is to multitype 
branching processes, and the dominant eigenvalue of the mean matrix of the 
branching process has a corresponding interpretation in the biological con- 
text. For more detailed discussions of such issues, see Heesterbeek (1992) 
and Diekmann & Heesterbeek (2000, Section 5.7). 

Whittle (1955) was able to justify his birth and death approximation to 
the early stages of the Markovian SIR-epidemic, and hence his formula for 
the probability of a large epidemic occurring, by sandwiching the epidemic 
process, during its initial stages, between two birth and death processes with 
slightly differing transition rates. This can be interpreted in terms of a 
pathwise comparison of processes. Ball (1983) and Ball and Donnelly (1995) 
went rather further, using a coupling argument to link the epidemic process 
with an approximating branching process on one and the same probability 
space, in such a way that the paths of the two processes are identical for a 
certain (random) length of time. In particular, they showed that the total 
variation distance between the distributions of the paths of the branching 
and epidemic processes is small, up to the time at which M = infections 
have taken place, for any choice = o(y/~N). They also suggest that this 
range of M N cannot be extended. 
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The coupling used by Ball and Donnelly is simple and natural, and it 
is somewhat surprising that accurate coupling is in fact possible, for some 
epidemic processes, over rather longer time intervals than they had supposed 
possible. This was first established by Barbour and Utev (2004), in the 
context of the discrete time Reed-Frost epidemic process. They showed that 
the branching process approximation to the path distribution actually has 
asymptotically small error in total variation for all choices of = o(iV 2 / 3 ). 
The essence of their argument lay in examining the likelihood ratio of the two 
processes along paths of given length, and showing that it was typically close 
to 1. In this paper, we show that similar arguments can also be applied to 
some continuous time models. We take as example the infinite dimensional 
BK-model, introduced in Barbour and Kafetzaki (1993) and subsequently 
generalized by Luchsinger (2002a,b), for describing the transmission of the 
parasitic disease schistosomiasis. 

2 The BK-model 

In the BK-model, N hosts are infected by parasites, with X^t) hosts hav- 
ing j parasites at time t, for j G Z + and t > 0. The process evolves as a 
Markov jump process X N in continuous time on the set X := G Z + , j > 
0) : J2j>o^i = w ith transition rates given by 

f -> f + e(j - 1) - e(j) at rate j/j,^, j > 1; 

f -> £ + e(j)-e(0) at rate A&X^/^K' j > 1, 

i>i 

for any £ G X, where e(j) denotes the unit vector in the j-th coordinate. 
The first of the transitions models the death of a parasite in one of the £j 
hosts currently carrying j parasites, the parasites being assumed to have in- 
dependent exponentially distributed lifetimes with mean The second 
transition models infection. Only currently uninfected hosts can be newly 
infected, and each makes contacts that could potentially lead to infection at 
rate A, the chance of such a contact being made with a host carrying I para- 
sites being (£i/N) (homogeneous mixing of hosts). If there is such a contact 
between an uninfected host and an /-host, then j parasites are established 
in the previously uninfected host with probability pij] in the BK-model, it 
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is supposed that pij = F[Ui = j], for U x := Y^i=\^i-> where the are inde- 
pendent and identically distributed random variables with mean 9 and finite 
variance, implying that each of the I parasites transmits on average 9 infec- 
tive stages to the newly infected host at an infectious contact, independently 
of the others. 

For a disease such as schistosomiasis, infection is actually indirect, and 
involves a host infecting suitable aquatic snails and these snails subsequently 
passing infection to other hosts. Thus the BK-model does not seem at first 
sight to be at all realistic. However, it can be thought of as an extreme case 
of a model incorporating features of the transmission process that were not 
present in many of the previous models: infection by parasites in groups, 
rather than singly, immunity in the definitive host (here, in the form of 
perfect concomitant immunity), explicit incorporation of the parasite burdens 
of individual hosts. 

The model that results is interesting for a number of reasons. The first is 
that, although it is rather complicated, it is still simple enough for some ana- 
lytic conclusions to be reached. For instance, it can be shown that the model 
has a 'law of large numbers' approximation for large N, in the form of the 
solution to an infinite system of differential equations, whose components 
approximate the proportions of hosts with different numbers of parasites. 
If 9 > e, this differential equation system has no (endemic) equilibrium solu- 
tion that yields a finite mean number of parasites per host. In practice, the 
distribution of parasites among hosts is observed to be extremely irregular, 
so that such behaviour is very encouraging: most earlier models have tacitly 
predicted Poisson-like distributions, which are far from realistic, and those 
that have tried to account for the over-dispersed distributions observed have 
imposed a specific form for the distribution without proposing any mech- 
anism that might generate it. Another feature is that, if 9 < e, there is 
exactly one equilibrium distribution of the differential equation system that 
has finite mean number of parasites per host, and that, in this equilibrium, 
the distribution of the number of parasites per host, conditional on the host 
being infected, depends only on the value of 9, and not on A or fi. 

For the purposes of this paper, it is the behaviour when few hosts are 
infected that is of primary relevance, with interest centering on questions such 
as the probability that the introduction of a single infected host can cause 
endemic infection to become established. These are the kinds of problem 



4 



that can be addressed by way of a branching process approximation. Here, 
we begin by proving an error bound for the approximation (Theorem 13.11) 
that is asymptotically valid in total variation for paths of length o(iV 2//3 ) 
transitions as iV — > oo. The branching process in turn yields a criticality 
theorem, which, to a close approximation, describes whether or not endemic 
equilibrium is possible in the BK-model. 

However, the approximating branching process — a Markov branching 
process with countably infinitely many types — itself displays unexpected 
critical behavour. If 9 < e, the branching process is super-critical, in the 
sense of having positive probability of growing indefinitely, if and only if 
X9 1 n > 1. The quantity \9 j \i has an immediate interpretation, being the 
lifetime average number of offspring of a single parasite, where offspring is 
interpreted in terms of parasites successfully passed on to other hosts, and 
is therefore precisely the biological quantity i?o, as seen from the parasites' 
viewpoint. Its appearance as the criticality parameter is therefore exactly 
what one would expect. However, if 9 > e, the criticality parameter is 
Aelog^//!, a fact that is much more difficult to interpret. 

Another feature of the model is that the mean number of parasites de- 
velops in time with exponential rate X9 — fi, whereas, if 9 > e, a super- 
critical process has a smaller exponential growth rate for the number of in- 
fected hosts. Thus, in such circumstances, the mean number of parasites 
per host increases ever faster. As a result, because deaths of parasites are 
counted as transitions, paths containing transitions may contain many 
fewer infections — roughly speaking, one may well have only ~ 
infections, for some a < 1. For such choices of the parameters, this makes 
the above theorem unsuitable for direct comparison with the results of Ball 
and Donnelly (1995). We therefore prove a second error bound in Theo- 
rem 13.21 which is expressed in terms of the asymptotics of Mjy. Its proof 
turns out to be a relatively simple adaptation of that of Theorem 13.11 We 
conclude with Theorem I4.2[ which establishes a rather stronger local state- 
ment, showing that the ratios of the likelihoods under the two models of paths 
containing at most M infections typically differ from 1 by more than order 

O ( (M 2 / 3 /N)^\og(N/M 2 / 3 ) ) with asymptotically negligible probability. 
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3 Total variation approximation 



The Markov branching process X := (Xj(-), j > 1) that approximates the 
BK- model is obtained from the process X N by ignoring the 0-component, 
taking the countable set X* : = {(£,■ G Z+, j > 1): X^>i£? < °°} as state 
space, and modifying the transition rates to 

f -»■ £ + e(J - 1) - e(J) at rate j/i^-, j > 2; 

£ -»• £ - e(l) at rate 

£ -»■ C + e(i) at rate Aj^&p y , j > 1, 

for £ E X*. These rates are identical with those for X N , except that, in 
the infection transition, the factor £o/N = 1 — S(£)/N is replaced by 1, 
where S(£) := Ylj>i£j- This represents the fact that, in the branching 
approximation, the total proportion of infected hosts is considered to be 
vanishingly small. Clearly, this should make little difference to individual 
transitions if S(£) <C N. The main result of this paper is to show that 
it makes little difference even for the distribution of whole path segments, 
considered as paths in X*, provided that the number of transitions m in the 
segment and the initial state £ (0) G X* are such that m + S(^) < iV 2/3 . 
We denote such a path by {(^ l \t^), < / < m}, where := 0, and we 
let T m denote the Borel a-algebra of events generated by these paths. To 
avoid trivial exceptions caused by paths that are absorbed in G X* never 
making further jumps, we suppose that both processes, when in state 0, make 
'jumps' to state at unit rate. 

Theorem 3.1 Suppose that ^ G X* , N > 2 and m are such that S m ^/m < 
N, where S m := m + S(^). Then, for any A G T m , we have 

\F[X N * G A] -F[X G A}\ < 8^^, 

where X N * denotes the process X N without the zero coordinate. 
Proof. For f G X* with 1 < S(f ) < N, write 

A N (0 := A(l-S(0/JV)£6(l-Pio), P N (t) = »J2& + AN (^ 

i>i j>i 

A(0 := A 6(1 -Pio), P(0 = • Mi). 

1>1 j>l 
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The quantities p N (£,) and respectively denote the overall jump rates of 
the processes X N and X in state £, A N (C,) and A(£) the overall infection 
rates; for ( = 06 X*, we set p N (0) = p(0) = A 7V (0) = A(0) = 1. Suppose 
that m + S(£^) < N. Then, for a path with m transitions starting in 
at time and then passing through the sequence of states 1 < I < m) 
at times < t {2) < • • • < t (m \ the likelihood ratio d¥ X N/dF x evaluated at 
such a path is just 

■= n{ex P {-(A-(^)-A(^))(^-^)} (i-^yy.i) 

where 

U| fl if S&W) = 5(^0) + i; 
1 otherwise. 

Hence we have 
with 

< := exp{-(A*(£«) - A(^ ))(t a+i) _ t (,) )} _ ^ 



and 

Note that each of these quantities is zero if £^ = G X*. 
The inequality 

< < S(^)/N (3.2) 

is immediate. Then, from the definitions of A N and A, it follows directly that 

implying that 

|A N (0/A(0 - 1| < (3.3) 
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Furthermore, if 

fj% : = |A w (^)/A(e (0 )-l|A(e (0 )(t (m) -t (/) ) < 1, (3.4) 
we also have 

K| < 2\fj%\ < 2{S(e i) )/N}e h (3.5) 



where 



ei := p^ l W l+1) -t {l) ) 



Hence, if (13. 4 p is satisfied, it follows that 

\L? +l -L?\ < N- 1 S(^){l + 2e l }L?. (3.6) 

Now suppose that (X®, I > 0) is a path resulting from a realization 
of the process X starting with X(0) = £ {0) , and that (T®, I > 1) are the 
corresponding jump times: set = 0. Then, defining 

E l+1 := p(X w )(T (i+1) -T w ), />0, (3.7) 

we note that £(Ei + i \ T\) is the standard exponential distribution for each /. 
Furthermore, the process {Lf(X^'\ "'), I > 0} is a non-negative martingale 
with Lq = 1 a.s., and so is the stopped version := N N (X('\ T^'>), 
where 

rf := inf{/>0: E l+1 > N/S(Z®)}, 
r-f := M{1 > 1: L? > 2}. 

Note also that S(X®) < S(£ (0) ) + / < S m for all < / < m. Now it 
follows from (13.31) and (13.71) that (13. 4p is satisfied for all < I < t± . Hence, 
from (13. 6p and the definition of t% , it follows that 

\L? +1 - Zf| < 2S m N- 1 (l + 2£ ; ) for aU < I < m. 

Thus, because Lf is a martingale, it follows that 

~ dm 9 2 ^9m 9 2 

E(L^-l) 2 < ^Ea + 2^) 2 = =^ (3.8) 

Now, for any A G jF m , we have 

P[X G A] - F[X N * G A] = E{(l-lf)l{Ievl}} < E{(l-0 + } 
< P[rf < m] + Ptfr-f < m} n {rf > m}] + E(l - Z*)+ (3.9) 
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From the definition of ly , it is immediate that 

P[rf <m] < me- N/s ™ < 4e" 2 (3.10) 

if S m ^Jrn < N. Then we have 

¥[{r 2 N <m}n {rf > m}} < - 1 > 1] < E(Z* - 1)+. (3.11) 

Finally, we have 

E(l-^) + + E(Z^-l)+ = E\l-L%\ 

< ^/e(1 -Z^) 2 < 2Vl3S m V^/N. 

It remains to note that 4e~ 2 + 2a/13 < 8. □ 

In general, the bound given in the theorem provides useful information as 
long as S m y/m N. In asymptotic terms, for fixed this allows paths of 
lengths mjv = o(iV 2 / 3 ) as iV — > oo, with an error bound of order 0(iV 37 / 2_1 ) 
if rriN ~ A^ 7 for some 7 < 2/3. 

If the Ball and Donnelly (1995) coupling is used to obtain error bounds, 
the resulting order 0(N 2 " , ~ 1 ), if is fixed and ~ A^ 7 , is at first sight 
not as sharp. However, there is an important difference between the two re- 
sults: the theorem above has m^, the total number of transitions, in the error 
bound, whereas the Ball and Donnelly coupling leads to an error expressed 
in terms of M^, the number of births or infections. Now the total num- 
ber of transitions includes all the parasite deaths, and if the mean number of 
parasites per host grows fast, as may be the case when 6 > e, may be sub- 
stantially bigger than M^. Thus Theorem 13.11 is not strong enough to yield 
an obvious improvement. For this reason, we now bound the discrepancies in 
the likelihood ratio more carefully, basing the argument explicitly on the se- 
quence of infection events. To this end, we let Hi denote the Borel a-algebra 
of events generated by paths containing exactly I infection events; as before, 
to avoid trivial exceptions caused by paths that are absorbed in G X* hav- 
ing no further infections, we suppose that both processes, when in state 0, 
create 'pseudoinfections' at unit rate. 



9 



Theorem 3.2 Suppose that^ e X* , N > 2 and M are such thatS M \fM < 
N , where Sm := M + S(^). Then, for any A E Hm, we have 



-N* ^ 41 TOtV r- /111 / qSmVM 



\F[X N * e A] -¥[X E A]\ < 8- 



N 



Proof. The likelihood ratio at a path £(•) in TCm can be written, using (13. ip . 
in the form 

L' N M := L%(Z(u), 0<u< a M ) 

■■= n{ex P f-£(A JV (ew)-A(ew))^ (i-^^)}(3.i2) 

where = a < o\ < • • ■ denote the times of infection transitions. Hence, 
very much as before, we have 

L l = il-lfl+lldKl-llB), 

with 

rfn ■= exp ( - P (A*(£(«)) - A(£(«))) du ] - 1, 



and 

Vl2 — ^ , 

each of these quantities being zero if f (0 = oer 
From (I3.3p . setting 



^ := f {i-A w (ew)/A(ew)}A(e(«))d« (3.13) 



we have 



< Vn < Wit < 2{S(£( i_ 1 ))/JV} C|) (3.14) 



~ AT 

whenever r/' n < 1, where 



ei := / A(f(u))du. 
10 



Hence, noting that S(£(u)) < Sm for < w < <jm, if 

ei = f A(£(u))du < N/S M (3.15) 

is satisfied, it follows that 

\Lf -L'ZA < N- 1 S M {l + 2e l }Lf_ 1 , 1 < I < M. (3.16) 

Now suppose that (X(u), u > 0) is a path resulting from a realization 
of the process X starting with X(0) = and that (<Ji, I > 1) are the 
corresponding times of births (infection transitions): set <r = 0. Then, 
defining 

E[ := f 1 A(X(u))du, I > 1, (3.17) 

we note that £(i? z / +1 | W;) is the standard exponential distribution for each /. 
We now argue as before using the likelihood ratio martingales {L'f r (X(-)), I > 
0} and L\ := Lf At ,N At ,N(X(-)), where 

t'f := inf{/>0: E' l+1 >N/S M }, 
:= inf{/ > 1: L'f > 2}. 

Since (13.151) is satisfied for all < I < t'f , it follows from (I3.16P and the 
definition of t'2 , that 

\L'" +1 - L'f| < 2S M N-\l + 2E' l+1 ) for all < I < M. (3.18) 

The remaining argument is now exactly as before, using the martingale L'i 
to compare the probabilities Pp^* G A] and F[X E A] for A G H A /- □ 

Thus Theorem 13.21 yields bounds of order 0(iV 37//2_1 ), improving on the 
rate obtained using the Ball and Donnelly (1995) coupling, if is fixed 
and Mtv ~ iV 7 for 7 < 3/2, where denotes the number of infection 
transitions. 

The new argument exploits the fact that the life histories of individuals 
infected with a given number of parasites have identical distributions in both 
models, except for the infection events, so that the likelihood ratio is corre- 
spondingly simpler. The key element is then that the difference in infection 
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rates between the two models is sufficiently small compared to the infection 
rate itself. The argument in Theorem 13.11 is less precise largely because, if 
the number of parasites is large, the bound (13.51) is rather pessimistic, since 
a potentially small factor A(£)/p(£) is not being exploited. 

4 Local approximation 

It was argued in Barbour and Utev (2004) that total variation approximation 
is not necessarily the best measure of closeness, if statistical applications in- 
volving likelihoods are to be justified. It is much more natural to want to have 
local approximations, which ensure that the ratio of actual and approximate 
likelihood is very close to 1, except possibly on a set of very small proba- 
bility. As a result, they defined a measure of relative closeness: probability 
measures P and Q on T are said to be e-relatively close with tolerance 77, 
RC (e, rj) for short, if there exists a set R G T such that 

P{R C ) < rj, Q{R C ) < 7], sup\log{{dP/dQ){x))\ < e. 

In this section, we show that the branching process approximation of the 
previous sections is indeed relatively close, as long as M N <C N 2 ^ 3 . 

We begin with a minor modification of the bounded differences lemma 
for martingales. 

Lemma 4.1 If(L n ,Q n , n > 0) is a martingale, and if 

\L n+ i — L n \ < a + bE n+1 for each n>0, 

where C(E n+ i | Q n ) is the standard exponential distribution exp(l) for each 
n > 1, then 

( 1) : max{P[L n - L < -y] ,F[L n -L >y}} < exp 
for all 

4:71 

< y < — e {(a + b) 2 + 6 2 }/ max(a, 6), 
3 

where eo > 1/15 is the constant defined by e eo (l — £o) -3 = 4/3. Furthermore, 
for all y > 0, 

(2) : max{P[L n -L < -y], ¥[L n -L > y}} < exp (- ~ V + -M . 

L 15 max(a, ojyn 135 J 



Sy 2 



ln{(a + b) 2 + b 2 } 
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Proof. If X is any random variable with EX = and \X\ < a + bE, where 
E ~ exp(l), then it follows that, for any 9 > 0, 

E{e ex } < E{l + ^ + imW} 

< l + ^ 2 E{(a + bE) 2 e^ a+bE ^} 

„ , n o n0 { a 2 2ab 2b 2 
1 + ^0 e ad < — + tz 7^ + 



2 

2 



i -be {i -bey (i-bey 



< exp<{^ 2 {(a + 6) 2 + 6 2 }^ 



as long as 6*max(a, b) < eo, with e defined as above. Hence, for any n > 1, 
E { e e(L "~ Lo) | £ n -i} = e e(L "- 1_io) E (e^^-^-i) | £ n _ x } 

< eX p |^ 2 {(a + &) 2 + & 2 }} e^- 1 "^, 

implying that 

E |^(l„-l )} < exp|V{(a + 6) 2 + 6 2 }|E{e e ( L "- 1 - L »)} 
for all n > 1, and hence that 

E { e 0(L„-L o )j < exp|^ 2 {(a + 6) 2 + 6 2 }|. 

Hence, for any y > and any e such that 6 max(a, b) < e , we have 

nL n -L >y] < exp|-^ + ^ 2 {(a + 6) 2 + 6 2 }|. 
Now, if |/max(a, b) < (4e /3)n{(a + b) 2 + 6 2 }, we can take 

a = V 

' ~ {4n/3){{a + b) 2 + b 2 } , 

to give 

p[i„-L„>yi < "p{ 8n{(a ; 3 g + ^} }- 
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On the other hand, for all y > and n > 1, we can choose 9 = 1/{15 max(a, b) y/n}, 
giving 

nL n -L >y] < exp(- ~ V + -M . 

[ 15 max(a, ojA/n 135 J 

The same arguments also cover P[L„ — Lq < —y] for the corresponding 
choices of y, since the conditions of the theorem apply equally well to the 
martingale — L n . □ 

This lemma enables us to prove the following estimate of relative close- 
ness. 

Theorem 4.2 Suppose that^ E X* , N > 2 and M are such thatij)(M, N) := 
SmVM/N < 1, where Sm '■= M+S^ ^). Then, with respect to paths in7i M , 
the processes X N * and X are RC (s r M N , Vmn) relatively close for any choice 
of r > I, where 

e r M,N ■= C r ip(M, N)y/log(l/ip(M, N)); 

rf MN ■= 2ip(M, N) r + e 2/135 exp{-l/60^(M, A^)} + Me~ N/SM , 
and C r := y/416r/3, provided that M > (l/5)C r 2 log and that e r MjN < 1. 

Proof. It was shown in the proof of Theorem 13.21 that the likelihoods of the 

processes X N * and X are close; here, we tighten the argument. We start 
from f)3.18p . which states that 

\L'" +1 -L'?\ < 2S M N- 1 (l + 2E' l+1 ) forallO</<M, 

where C(E' l+1 \ Hi) is the standard exponential distribution for each I, and 
from the observation that, by the definition of L' M , we have L' M = L M as 
long as minjt'f, t'^} > M. Now it is immediate, as for (I3.10p . that 

P[t'f <M}< Me~ N/SM . 

Then, from the definition of t'2 , it follows that 

P[{t'f > M} n {tf < M}] < ¥[D N M - 1 > 1]. 

Hence, to establish the desired relative closeness, we take 

R c := {min(tf ,tf ) < M} U {\L' M - 1| > e r M)N /2}, 
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(here using the assumption that e r M N < 1) and bound the probabilities 

P[L' M - 1 > 1] and F[\L' M - 1| > e r M , N /2] using Lemma O with n = M and 
2a = 6 = ASm/N. 

First, we use Lemma [4.11 (2) to give 

v „ n . r rrrs, 2/135 



F[L' M - 1 > 1] < exp{-iV/(605 A/ v / M)}e 
Then we use Lemma [4.11 (1) to show that 

¥[\1-L'm\>v] < 2exp{-3iVV/(416M^ / )} 



provided that 

4M 135 M 

< y < e — — — . 

_ y - 3 o N 

Hence we can take y = e r M N /2 if 



C r S M VMy/logN/N < VQAMS M /45N 
and thus if M > (l/5)C^logiV, giving 

. >v (3/416) ( 



l'1\> --:;,.v/2] < - ; ^ > = 2<j^pi> . □ 



Thus asymptotic relative closeness of order 0{ifj(M } N) ^log(l/tp(M, N))} 
can be established with tolerance of arbitrarily small polynomial order in 
ip(M, N) = SmVM/N. 
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